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ABSTRACT 

/ 

'\ 

We  consider  the  choice  of  technique  in  a  continuous  time  infinite 
horizon  optimal  growth  model.  There  are  n  +  2  goods,  output,  labor 

and  machines  . . . ,  M'  .  We  can  convert  one  unit  of  labor  to 

r 

q  units  of  output  or  r^  units  of  ,  for  each  i  ,  Also,  we  can 

convert  one  unit  of  labor  and  one  unit  of  to  q^  units  of  output. 

We  prove  under  some  sufficient  and  necessary  conditions  that  we  never 
build  any  machines  for  the  general  concave  utility  function.  If  the 
condition  is  not  met,  we  build  one  machine  from  beginning  to  end  when 
the  utility  function  is  linear;  when  the  utility  function  is  non¬ 
linear  life  gets  complicated.  In  the  one  machine  case,  we  give  a 
general  algorithm  to  solve  it.  In  the  many  machines  case,  we  prove 
an  asymptotic  result  (as  t  ^  ,  the  behavior  is  similar  to  that  of 

the  linear  case)  and  give  examples  showing  that  a  simple  characteriza¬ 
tion  of  the  optimal  solution  is  difficult.  , 
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CHAPTER  I 


MODEL  AND  DISCUSSION 


II.  Introduction 


In  this  paper  we  investigate  the  choice  of  investment  in  a  con¬ 
tinuous  time  optimal  growth  model.  The  general  model  is  described 
as  follows:  there  are  m  goods,  some  of  which  may  be  provided 
exogenously  (e.g.  labor).  The  technology  is  described  by  a  set  of  n 
activities,  each  of  which  consumes  various  amounts  of  goods  and  produces 
various  amounts  of  goods.  To  be  more  specific,  the  technology  is 
given  by  a  pair  of  nonnegative  m  x  n  matrices  A  B  and  a  nonnegative 
n-vector  b  ,  where  denotes  the  amount  of  good  i  used 

(produced)  to  operate  activity  j  at  unit  level,  and  b^(t)  denotes 
the  amount  of  good  i  exogenously  provided  at  time  t  .  There  is  a 
utility  function  which  is  an  increasing  concave  function  of  the  activity 
level. 


Problem: 

Given  b(t)  ,  find  an  activity  vector  x(t)  to  maximize  the  dis¬ 
counted  integral  of  future  utility. 

We  can  write  this  as  a  continuous  programming  problem: 


<P> 


Maximize 


J  e  °''*'*U(x(t))dt 

*0 


subject  to:  A*x(t)  <  I  B*x(s)ds  +  b(t) 


1' 


x(t)  >  0  . 


■"•r  I  lift'll  iTi  iTi'iiTi 
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If  x(t)  is  <P>  feasible  ^x(t)  >  0  ,  Ax(t)  <  J  Bx(s)ds  +  b(t)| 
then  we  have  the  following  dual  problem: 


<D(x(t))>:  Minimize 


00 

w(t)  • 


b(t)dt 


00 

/“( 


subject  to:  w(t)A  >  J  w(s)Bds  +  e  *^^*VU(x(t)) 

t 


w(t)  >  0  . 


§2.  Optimality  Theorem  for  the  General  Model 

In  this  section  we  prove  that  a  feasible  solution  (x(t))  is 
optimal  if  we  can  find  the  corresponding  dual  price  (w(t))  and  satisfy 
some  complementary  slackness  conditions. 


Theorem  1-1: 

If  the  following  four  conditions  are  satisfied 


(1) 

(2) 

(3) 

(4) 


x(t)  is  <P>  feasible  /i.e.,  Ax(t)  <  /  B*x(s)ds  +  b(t)  , 

i(t)  >  ol  . 

D(x(t))>  feasible  li.e.,  w(t)A  >  J  w(s)Bds  + 

\  V  “  t 


w(t)  is  < 
-at 


VU(x(t))  ,  w(t)  >  0 


) 


w(t)^Ax(t)  -/  B*x(s)ds  -  b(t)J  =  0 
|^(t)A  -  J*  w 


(s)Bds-e  ‘VUCxCt))  ‘xCt) 


I  ' 


then  x(t)  is  <P>  optimal,  w(t)  is  (D(x(t))>  optimal. 
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Proof; 


Let  x(t)  be  any  (P)  feasible  solution 

w(t)  be  any  <D(x(t)))  feasible  solution. 


w(t)  'b ( 


t)  >  w(t)*  A*x(t)  -J  B-x 


(s)ds 


Integrating  both  sides  gives 


w(t)*b(t)dt  >  J  w(t)-|^Ax(t)  -  J*B*x(s)ds|dt 

00  "  OO  “ 

=  I  w(t)A  -  I  w(s)Bds  • 

*0  L  *c 


w(s)Bds  •x(t)dt  (changing  the  order  of 
integral) 


>  I  e  •VU(A(t))*x(t)dt  (w(t)  is  <D>  feasible) 


>  I  e  [VU(x(t))*x(t)  +  U(x(t))  -  U(x(t))]dt 


00  00 

J  e~°‘^u(x(t))dt  +  J'  e‘°‘''[VU(x(t))-x(t)  -  U(x(t))]dt 


e  •U(x(t))dt  +  M  . 


So  for  any  x(t)  ,  w(t) 


w(t)-b(t)dt  >  I  e~“’^*U(x(t))dt  +  M 


but  if  we  replace  x(t)  ,  w(t)  by  x(c)  ,  w(t)  then  the  above  ”>" 
are  replaced  by  (from  (3)  and  (4))  thus 


00 


1 


w(t)*b(t)dt  “  I  ®  '*^*U(x(t))dt  +  M 

'0 


From  (5)  and  (6) 


00 

I 

00 


w(t)-b(t)dt  reaches  the  lower  bound  and 


e  ^^•U(x(t))dt  reaches  the  upper  bound. 


So  w(t)  is  <P)  optimal 

w(t)  is  <D(x(t)))  optimal# 

§3.  Description  of  the  Investment  Model 

We  now  consider  the  special  case  as  follows; 


(1)  There  are  n  +  2  goods,  output  (consumption),  labor,  and 

n  machines  M  ,M  ,  ...,  M 
1  Z  n 

(2)  There  are  n  +  1  production  activities  P  ,P  ,  ...,  P 

0  1  n 

Pq  converts  one  unit  of  labor  to  q  units  of  output. 

P .  converts  one  unit  of  labor  and  one  unit  of  to  q . 

i  i  1 

units  of  output  1  <  i  <  n  . 

(3)  There  are  n  investment  activities  1,  ,1-,  ...,  I  ,  where 

12  n 

converts  one  unit  of  labor  to  r^  units  of  ^  i  j 

(4)  Labor  is  given  exogenously  at  constant  rate  of  one. 


Let  X  be  the  amount  of  labor  allocated  to  ; 


z  be  the  amount  of  labor  allocated  to  P .  1  <  i  <  n  ; 

i  1  -  - 

y  be  the  amount  of  labor  allocated  to  I.  1  i  i  i  • 

^  i  1  -  - 

Let  a  be  the  discount  rate  and  U  be  a  differentiable  increasing 
utility  function  of  the  output. 

n 

We  define  c(t)  =  q*x(t)  +  [  q  -z  (t)  and  tp(t)  =  U'(c(t))  . 

i=l  ^ 


Problem; 


Given  initial  stocks  of  ,  choose  X(t)  -  (x(t),z^(t), 


y^(t))  e  to 


! 


-at 


maximize  I  e  •U(c(t))dt 

*0 


subject  to:  z 


i<'>  i  J >'i 

^r\ 


(s)ds  1  1  ^  5  ^ 


n 


x(t)  +  \  (y.(t)  +  z  (t))  <  1 

i=l  ^ 


By  defining  U(X(t))  =  U(c(t))  ,  VU(X(t))  =  (q,q^,  •  •  • ,  q^ »  0>  •  •  •  »  0)  *  ^  (t^) 

n  terms 


^  TD^n+l  j 

G  ]R_^  and 


0  ;  1 

0 

1  ’'l 

• 

1  r_  0 

s 

■  1  I  1  0 

,  B  = 

• 

0  1 

0  ‘  * 

^  1  _1 

0 

1  0  ‘^n 

ill,  ...«  ill,  •••»  1 

0 

0 . 0  I  0,  . . .  ,  0 

A  .  B  e  ]R("-*-l>"(2n+l)  ^  k  .1)  ,  b(t)  G  TR] 


n+1 
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The  dual  problem  becomes: 


n+1 


Find  w(t)  =  fw  (t),  ....  w  (t),w(t))  G  IR  to 
in  + 


00 

1 


minimize  I  w(t)*b(t)dc 

*0 

"“Ctt 

subject  to:  w(t)  >  e  -q'^t/Ct) 


d^(t)  =  w(t)  -  r.  I  w. (s)ds  >  0 


00 

r.  I  w.  ( 


1  <  i  <  n 


-»ott 

w(t)  +  w^(t)  >  e  'q^'il^Ct) 


1  <  i  <  n 


If  the  utility  function  is  linear,  by  giving  the  exact  dual  price, 

,  the  stationary 


i  ^ 

we  prove  that  (1)  if  q  >  max  \ — '•  }  =  — ; — 

^=,^.^la  +  r.l  cx  +  r 

program  (y^(t)  =  0  ,  1  <  i  <  n)  is  optimal;  (2)  if  q  <  ^  , 

/  \  ^ 

we  build  M  (in  fact,  any  M  s.t.  — - - -  — - - )  .  If  the 

Z  [  ’  ^  t  a  +  r^  a  +  r J 

utility  function  is  nonlinear,  life  gets  complicated. 

To  begin  with,  we  do  not  know  whether  a  right-differentiable 

(see  Perold  [5])  optimal  solution  exists.^  Secondly,  we  do  not  know 

2 

whether  the  dual  price  exists.  To  simplify  the  model,  we  restrict 
the  utility  functions  to  those  "well-behaved"  differentiable  increasing 
functions  such  that  there  exist  right-differentiable  optimal  solutions 


Since  our  model  is  an  economic  model,  we  would  like  an  optimal  solution 
which  is  right-differentiable  rather  than  just  Lebesque-measurable. 

In  1968,  Hanson  [3]  proved  the  strong  duality  of  the  similar  model  in 
the  finite  horizon  case. 
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and  the  corresponding  dual  prices.  Under  this  assumption,  we  prove 

that  (1)  q  >  ^  ^  is  still  the  sufficient  and  necessary  conditions 

for  the  existence  of  an  optimal  stationary  program;  (2)  if  the  condition 
of  (1)  is  not  met,  we  can  prove  an  asymptotic  result,  i.e.,  as  t  ®  , 
the  optimal  solution  has  the  same  behavior  as  the  linear  case. 

It  is  difficult  to  give  a  simple  characterization  of  the  optimal 
solution  of  a  nonlinear  utility  function  because  we  have  examples  in 
which  we  switch  building  machines  (Chapter  4,  Example  1),  or  build 
some  machine  in  time  [1^,12]  but  later  never  use  . 

Even  in  the  one  machine  case,  we  might  have  x(t)  =  0  in  the 
first  interval,  x(t)  >0  in  the  second  interval,  and  x(t)  =  0 
in  the  third  interval  (Chapter  3,  Example  2). 

In  the  following  we  prove  that  by  manipulating  the  model  we  can 
make  some  assumptions  without  losing  generality. 


Assumption  1:  q,  >  q„  >  . . .  >  q  ,  r,  <  r„  <  . . .  <  r 
- -  12  n  1  2  n 

If  “  'Ij  snd  ’  apparently  we  never  build  (for 

has  the  same  output,  but  larger  • 


(i) 

(ii) 


If  k.  =  0  ,  then  we  discard  M.  . 

J  3 

If  kj  >  0  ,  we  define  a  new  problem  with  n  -  1  machines 

(without  M.)  and  the  initial  stocks  k  given  by 
J  ® 


k.“k.  +k.  ,k  =k 
1  1  j  s  s 


l<s<n,s?ti,j  . 

It's  easy  to  see  that  the  new  problem  is  in  fact  the  same 
as  the  old  one. 


3 

It  is  believed  that  for  all  "usual"  concave  functions,  these  properties 
are  true.  Since  our  model  is  a  simple  economic  model,  there  should 
exist  the  dual  price  and  a  well-behaved  optimal  solution. 
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■* 


Hence,  without  losing  generality  we  can  assume  ^  ‘Ij  » 


if  i  ^  j  . 


Now  we  rearrange  the  machines  in  order  of  decreasing  output. 


Assumption  2:  K  =  ^  k.  <  1 

i=l  ^ 


n  i  j+1 

If  ^  k  >  1  ,  let  ^  k  <  1  ,  ^  k  >  1  . 

■_1^  1-^  j_i-^ 


We  will  show 


that  the  general  problem  is  in  fact  equivalent  to  this  special  case. 


(1)  We  never  build  M  ,  s>j+l  /  j+1 

®  (for  y  k  >  1  and  q,  >  q-  >  . .  .  >  q  | 

z,a)  =  o  ,  ,>3^2  vtjo  .  \  lii  12  „j 

(2)  By  (1)  we  discard  ,  s  >  j  +  2  . 

(3)  Consider  a  new  problem  with  machines  and  q 

replaced  by  • 

By  (2)  and  (3),  it  is  clear  that  the  new  problem  is  essentially  the  same 

as  the  old  one,  but  |  k.  <  1  for  the  new  problem. 

i=l 

Before  ending  this  chapter,  we  prove  two  theorems  which  play  the 
key  role  in  the  following  chapters. 


Theorem  1-2: 


(1)  x(t)  +  y  (y . (t)  +  z.(t))  =  1  for  all  t  (full  employment). 

i=l  ^ 

(2)  There  always  exists  the  dual  price  w^(t)  =  |e  ^'^•q^*ti'(t)  -w(t)j^ 
1  <  i  <  n 


which  satisfies  the  conditions  of  Theorem  1-1. 
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Proof : 

n 

(1)  If  x(t)  +  I  (y.(t)  +  z  (t))  =  P(t)  <  1 

1*1 

since  x(t)  has  output  q  ,  we  can  increase  x(t) 
by  1  -  p(t)  to  produce  more  output,  contradiction. 

(2)  w^(t)  >  0  and  w(t)  +  w^(t)  >  e  “'■•q^*i|)(t) 
implies  w^(t)  >  |e  q^*!i)(t)  -  w(t)J^  . 

Define: 


k^(t) 


■  •  -I 


yj^(s)ds 


If  >  0  ,  Case  i:  2^(t)  <  k^(t) 

by  complementary  slackness  w^(t)  *  0  = 

[e  “'^•q^*ij<(t)  -  . 

Case  ii:  z.(t)  =  k. (t) 

1  i 

w(t)  +  w^(t)  =  e  '^^•q^’il'(t)  implies 
w^(t)  =  je  “’^•q^*ii^(t)  -  w(t)J^  . 

If  k^  *  0  ,  Case  i:  there  exists  T  >  0  ,  such  that  k^(t)  ■  0 

for  t  <  T  and  k^(t)  >  0  for  t  >  T  • 
w^(t)  =  ^e  “’^•q^*ij'(t)  -  Vt  >  T  . 

Define  w^(t)  *  ^e  -  w(t)J^  for  t<T  . 

00 

w(t)  ^  f  Wj^(s)ds  (w^(t)  is  the  dual  price) 

for  t  <  T 

00  A  A 

>  rf  /  w^(s)ds  (w^(s)  >  W^(s))  . 


j 

i 
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y^(t)  =  0  for  t  <  T  (k^(t)  =  0)  implies 
the  complementary  slackness  is  also  satisfied 
for  w^(t)  . 

Case  ii:  if  k^(t)  =  0  for  all  t 
similar  proof  as  above. 

Thus  there  exists  an  optimal  price  w^(t)  *  fe  -  w(t)l_|_B 


Remark : 

(1)  If  k^  =  0  ,  w^(t)  need  not  be  unique. 

[— cic 

e  -  w(t) 


A  solution  X(t)  =  (x(t)  ,y^(t)  ,2^(t))  €  ]R_|_  is  called  stationary 


if 

x(t)  = 

1-k  ,  y^(t)  »  0  ,  z^(c)  “  k^  1 

k  »  I  k 

1  i=l  ^ 

j  for  all  t 

1  < 

i  <  n  , 

,  Let  max  — : -  »  v  ♦ 

^  V* 

l<i<n  1 

Theorem  1-3: 

q  >  V  is  the  sufficient  and  necessary  condition  for  the  existence 
of  an  optimal  stationary  program. 


Proof ; 

If  <1  ^  V,  the  stationary  program  is  feasible  and  c(t)  = 
n 

q*(l-k)  +  'l  k  'q  =  constant  implies  lii(t)  =  constant  =  c  >  0  . 
i-1  ^ 

Define:  w(t)  =  e  ‘q’lpCt)  “  e  'q'C  and  w^(t)  =  e  •  (q^  -  q)<('(t)  = 

-at  ,  . 

e  •(q^-q)*c  . 

It  is  easy  to  see  that  if  we  can  prove  d^(t)  >0  I  ^  I  ^ 
then  the  stationary  program  is  optimal. 


d^(t)  = 


> 


-at 

e  •q'c 


-  r. 
1 


00 


-as 

e  * 


q) 'Cds 


-at^ 

^  f“q  -  r^(q^-q)] 


aq 


r.(q.-q) 


is  equivalent  to  q  > 


“  '1/ 


The  other  direction  is  similar.  (Here  we  use  the  assumption  of  the 
existence  of  the  dual  price.)® 
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CHAPTER  II 


THE  LINEAR  UTILITY  FUNCTION 


To  develop  better  understanding  of  the  general  concave  utility 
function,  we  start  with  the  linear  utility  function.  There  are  only 
two  kinds  of  the  optimal  program  for  the  linear  utility  function: 
one  is  never  building  any  machine,  the  other  is  building  only  one 
machine. 

For  completeness,  we  list  the  stationary  program  in  the  following. 


§1.  Stationary  Program 


Theorem  2-1: 

q  >  V  is  the  sufficient  and  necessary  condition  for  the  existence 


of  an  optimal  stationary  program. 
From  now  on,  assume  q  <  v  . 


§2.  Nonstatlonary  Program 


Machine  j  is  called  best,  if  — +r 

N  ^ 

1<  i<n}  =  (1,2,  ...,  N}  ,  .  Machine  i  is  called  good, 

if  1  G  I  . 

For  any  best  machine  j  ,  the  following  program  is  called  the 
best  program  with  respect  to  j  :  x(t)  =  0  ,  y^(t)  *  0  i  j  , 

y^Ct)  -  (l-k)e  J  ,  z^(t)  -  k.  i  e  1  -  {j}  ,  Zj(t)  *  + 


(1-k) 


|l-e  ^  )  »  2^(t)  *  0  r  i  I  .  (This  means  that  all  the  good 


machines  are  fully  utilized,  then  the  rest  of  the  labor  is  used  to 


build  machine  j  .) 
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Theorem  2-2: 


The  best  program  with  respect  to  j  is  optimal. 


Proof : 


c(t)  =■  'I  q.‘z  (t)  ,  iKt)  =  constant  =  a  >  0  (because  U  is  linear), 
i=l  ^ 


Define: 


w(t)  =  e  '^^•[q*a  +  s]  =  e  ^  a  ,  S 

a  +  r . 


-  ■]■ 


a  >  0 


w.  (t)  =  e  • 

I  1 


w^(t)  =  e 


dj  (c)  =  e 


(q. - I  'a  if  i  €  I 

0  if  iCi  . 

f  -as  r  1 

•a*  — T"'* - I  s  *3'  hj,r - 1 

a^r.  jj  a  +  rj 


•a  =  0  if  i  C  I  . 


ds  =  0 


OD 

i  e  I-(j}  ,d(t)  =  e  °‘*'»a*  — - r .  T  e  '^^•a*  /q . - j;  \ 

i  a  +  r^  X  J  a  +  j 


-at 

e  ’a* 


— ii  >  0 

a  |a  +  r.  a  +  r.l* 

L  J  ij 


i  ^  I  ,  d^(t)  “  w(t)  -  0  >  0  . 


If  i  G  I  ,  z^(t)  -  k^(t)  ;  if  i  ^  1  ,  z^(t)  »  0  . 


The  complementary  slackness  is  also  satisfied 
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CHAPTER  III 

GENERAL  CONCAVE  UTILITY  FUNCTION  -  ONE  MACHINE  CASE 


For  the  general  concave  utility,  to  get  deeper  insight  about 
the  structure  of  the  many  machines,  we  begin  with  one  machine  case. 

In  the  beginning,  we  derive  some  simple  properties.  Later  we  prove 
that  if  the  initial  stocks  are  sufficiently  large,  we  don't  use  labor 
alone  to  produce  output.  If  the  utility  function  is  quadratic,  we 
have  a  simple  optimal  solution.  Otherwise,  we  have  an  example  with 
a  non-simple  solution.  We  also  give  an  algorithm  to  solve  the  general 
problem.  By  scaling,  we  can  assume  q  =  1  without  losing  generality. 

Because  of  the  simple  structure,  we  rewrite  the  primal  and  dual 
feasibilities  in  terms  of  the  following  simplified  notation: 


x(t)  x(t';  ,  y(t)  y^(t)  ,  z(t)  z^(t)  ,  r  *  r^  ,  q  q^^ 


z(t)  5  r  I  y(s)ds 


primal  feasibility: 


I 


x(t)  +  y(t)  +  z(t)  =  1 


x(t)  ,  y(t)  ,  z(t)  >  0 


at 


w(t)  >  e  •li'Ct) 


dual  feasibility: 


00 

/”i' 


'  w(t)  >  r  I  w  (s)ds 


w(t)  +  w^(t)  >  e  '^^•q'tpCt) 


w(t)  ,  Wj^(t)  >  0  . 
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For  completeness,  we  list  the  stationary  program  in  the  following. 

§1.  Stationary  Program 
Theorem  3-1: 

1  >  — —  is  the  sufficient  and  necessary  condition  for  the 
=  a  +  r 

existence  of  an  optimal  stationary  program® 

From  now  on,  assume  1  <  — —  . 

a  +  r 

§2.  Nonstationary  Program 

We  begin  with  some  simple  properties. 

Property  1: 

y(t)  >  0  ,  for  all  t^  . 

Proof : 

If  there  exists  an  interval  (T,T  +  e)  such  that  y(t)  =  0  for 

t  €  (T,T  +  e)  ,  then  y(t)  =  0  for  all  t  >  T  ,  which  by  stationary 

property  implies  1  >  — —  .  Contradiction  to  1  <  — ^ —  B 
f  j  r  =  a  +  r  a  +  r 


Corollary : 

00 

w(t)  -  r  /  w^(s)ds  =  0  for  all  t  B 


We  disregard  those  isolated  time  points.  When  we  say  y{t)  >  0  for 
all  t  ,  we  mean  there  is  no  interval  such  that  y(t)  =  0  on  that 
whole  interval. 
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Property  2: 


k.(t)  =  kg  +  r  /  y(s)ds  =  z(t)  (fully  utilize  the  machine), 


Proof : 


By  contradiction.  If  not  true. 


Case  i:  there  exists  T  >  0  such  that  z(t)  <  k(t)  for  all  t>T 

w^(t)  =  0  (complementary  slackness)  implies  w(t)  =  0  . 
w(t)  >  e  °‘^•t^l(t)  implies  i|)(t)  =  0  ,  contradiction. 


Case  ii:  there  exists  ,  T2  such  that  z(t)  <  k(t)  for 

t  G  (T^.T^)  and  z(t)  =  k(t)  for  t  £  (T2,T2  +  e;)  . 

C(T^)  =  x(T^)  +q-2(T^)  =  0+q.z(T^)  <  x(T2)  +  q*z(T2)  = 
0(12)  implies  t(T^)  >  iJj(T2)  . 

CO  00 

0  =  wd., )  -  r  /  w.  (s)ds  “  w(T. )  -  r  /  w  (s)ds  (w,.  (t)  »  0 
1  T  T 

^  ^  t  e  (T  ,12)) 

=  w(T^)  -  w(T2)  . 

-«Ti  _ 

0  <  w(T^)  -  e  ^-q-iKT^)  .  (w^(T^)  =  0) 

-«T2  _ 

w(T2)  <  e  ^’q‘<p(T2)  .  (zCT^)  >  0) 

"“'^2  -  "'^^l  - 

From  the  above  three  equations,  we  have  e  •q’'J'(T2)  ^  e  ‘q’lj^dj^)  •, 

contradiction  to  <  T2  and  '^(Tj^)  >  il'(T2)  B 

Corollary : 

k(t)  <  1  for  all  t  .  (We  never  accumulate  to  1  unit  of  machine.) 


Proof ; 


The  maximum  amount  of  machine  we  can  have  in  an  optimal  solution 


is  given  by  the  following  program:  x(t)  “  0  ,  y(t)  *  l-k(t)  , 
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z(t)  =  k(t)  .  k(t)  =  r  J*  [l-k(s)]ds  +  k^  ,  which  implies  k(t) 
1  -  (l-kQ)e  ,  thus  k(t)  <  1  for  all  t  ■ 


Remark : 


Property  2  does  not  hold  in  the  many  machines  case. 


Property  3: 


w(t)  is  continuous. 


Proof : 


w(t)  =  r  /  w^(s)ds  ■ 
t 

Property  4: 

li'Ct)  is  continuous. 


Proof : 


By  Property  2  and  i>(t)  =  U'(x(t)  +  q*z(t))  ,  if  we  can  prove 

for  any  given  T  >  0  ,  a  sequence  t^  -*■  T  such  that  * 

we  have  ijjCt  )  -►  i(j(T)  ,  then  <p(t)  is  continuous, 
n 

Case  i:  x(T)  *  0  .  If  x(t  )  ->■  x  =  0  ,  apparently  i|>(t  )  -*■  ’.>(T) 

n  n 

If  x(t  )  -*■  X  >  0  ,  then  lim  c(t  )  >  c(T)  and 
n  n  = 

n-Hx> 

limi()(t)<iJ'(T)  . 
n  = 

n-K» 

By  slackness,  e  )  =  w(t  )  -*■  w(T)  >  e  . 

n  n  “ 

Combining  the  above  two  equations,  we  have 

Case  ii:  x(T)  >  0  .  If  ^  same  proof  as  above. 

If  x(t  )  *  0  for  all  n  ,  lim  c(t  )  <  c(T)  and 
n  n  “ 

n-Hx> 


4 
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lim  >  iJ^(T)  . 

n-+«o 

e  <  w(cn)  — *■  w(T)  =  e  . 

Combining  the  above  two  equation,  we  have  "*■  '^(T)B 

Corollary  1: 

If  u(x)  is  strictly  concave,  then  (1)  c(t)  is  continuous; 

(2)  x(t)  ,  y(t)  are  continuous. 

Proof : 

(1)  If  c(t)  is  not  continuous,  then  is  not  continuous, 

contradiction. 

(2)  c(t)  =  x(t)  +  q"2(t)  .  For  c(t)  and  2(t)  continuous, 
x(t)  is  continuous.  y(c)  =  1  -  x(t)  -  z(t)  implies  y(t) 
is  continuous  tooB 

Corollary  2: 

Wj^(t)  is  continuous. 

Proof : 

w^(t)  “  ■•q'4'(t)  -  .  Since  li'(t)  and  w(t)  are 

both  continuous,  Wj^(t)  is  continuous  tooB 

Property  5: 

If  x(t)  >  0  ,  y(t)  >  0  ,  z(t)  >  0  ,  then  iJiCt)  =  c'e^*^  DH 
c  >  0  ,  t  e  (T^,T2)  . 


Proof : 


The  dual  equations  are;  w(t)  =  e  ,  w^(t)  =  e  -(q-Dii/Ct)  , 

w(t)  =  r/  Wj^(s)ds  .  The  solution  is  ip(t)  =  c-e^^  ,  c  >  OB 

t 

Corollary  1: 

If  u'(q)  >  0  and  x(t)  >  0  ,  y(t)  >  0  ,  z(t)  >  0  ,  then 
4j(t)  <  6  <  0  . 

Proof : 

ij(t)  =  [a  -  r(q  -  1)  ]'4/(t)  <  U  ’  (q)  [a  -  r  (q  -  1)  ]  =  B  <  OB 
Corollary  2: 

If  ij;(t)  =  constant  for  t  €  (T.T  +  t)  ,  then  x(t)  =  0  for 
t  G  (T,T  +  £)  B 

Corollary  3: 

c(t)  is  increasing  in  (T^,T2)  • 

Proof : 

Mt)  =  c-[a-r(i-l)]-e^“'’'^'^"^^^^  <  0  ,  so  ^(t)  is  decreasing  in 
t  .  Since  U' (x)  is  decreasing  in  x,  c(t)  is  increasing  in  (T^,T2)B 

Property  6: 

If  the  utility  function  is  strictly  concave,  then  c(t)  is  a 
strictly  increasing  function. 

Proof : 

If  there  exist  T^^,  T2  <  T2  such  that  c(Tj^)  >  c(T2)  , 

by  Property  A  there  must  exist  T  ,  T  G  [Tj^,T2)  such  that 
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c(T)  =  max  c(t)  . 

tG[T^,T2] 

If  x(T)  >  0  ,  by  Property  5  c(t)  is  strictly  increasing  in 

(T-e,T  +  £)  ,  for  e  small  enough,  contradiction  to 

c(T)  =  max  c(t)  . 

t€[T^,T2] 

If  x(T)  =  0  ,  by  Property  2  c(t)  is  strictly  increasing  in 
[TjT  +  e)  ,  for  e  small  enough,  contradiction  ■ 

In  the  following,  we  study  the  asymptotic  behavior. 

Theorem  3-2: 

If  the  utility  function  is  nonsaturate  (U'(q)  >  0)  ,  then  for 
t  sufficiently  large,  x(t)  =  0  . 

Before  proving  this  theorem,  we  need  the  following  lemmas. 

Lemma  1 : 

If  the  utility  function  is  nonsaturate,  is  the  initial  stock 

and  u'(q*kQ)  <  U' (q)  ♦  ^  .  Then  x(t)  =  0  ,  y(t)  =  (l-kQ)e  , 

z(t)  =  1  -  (1  -  kQ)e  is  optimal. 

Proof : 

Def ine 

(w(t)  =  e  “^•[i);(t)  +  s(t)] 

|w^(t)  =  e  (q  -  l)i|)(t)  -  s(t)] 

00 

where  s(t)  =  +  rq  J*  e  , 

t 
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00 

/ 


00 

>  (v(t)  is  decreasing) 


-ctt  -  ..s  a  +  r-  r.„ 

=  e  •q'lJ^Ct)  - - : - >  0  . 


a  +  r 


Apparently,  w(c)  >  0  too. 

oo 

By  calculation,  w(t)  =  r/w^  (s)ds  for  all  t  ,  which  implies 

t 

the  complementary  slackness  is  satisfied. 

If  we  can  show  w(t)  >  e'^'^'iKt)  ,  then  the  dual  feasibility  is 

also  satisfied,  thus  optimal. 

u'  =  u'  (q) 

00 

,  s  -  r  -(a+r)x,,  (r+a)t 

s(t)  =  ~ip(t)  +  rq  I  e  '  ip(x)dx*e 

t 

>  -U’(i-kQ)  +  7^-  U'(i)  >  0 


implies 


w(t)  >  e  *^*'’iJ;(t)B 


Lemma  2 : 

k(t)  -►  1  for  t  sufficiently  large. 


Proof ; 


Since  k(t)  is  increasing  and  bounded,  k(t)  -*■  K  ,  y(t)  -*■  0 


as  t  “ 


> 
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If  k  <  1  ,  as  t  ,  x(t)  =  l-y(t)-z(t)  =  l-k(t)-y(t)  -► 

1-k  >  0  .  By  Property  5  ij/(t)  -►  0  ,  contradiction.  Thus  k(t)  -»■  1 
for  t  sufficiently  large  ■ 

Now  we  are  ready  to  prove  the  theorem. 

Proof : 

^  ^  >  1  and  U'  (q)  >  0  imply  the  existence  of  k  <  1  such  that 
U’  (i-k)  <  U'  (i)  • 

By  Lemma  2,  there  exists  T  >  0  such  that  k(T)  >  k  . 

Now  by  Lemma  1,  for  t  sufficiently  large,  we  have  x(t)  =  OH 

In  the  next  section,  we  give  a  general  algorithm  to  solve  the 
one  machine  problem  completely.  This  algorithm  comprises  of  solving 
two  subproblems  (one  with  three  activities,  the  other  with  two 
activities)  backward  by  turns.  It  turns  out  that  we  can  partition  the 
machine  stock  into  disjoint  intervals  such  that  on  each  interval 
either  two  activities  or  three  activities  are  optimal,  and  alternately. 
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PROBLEM  I  (k) 

(1)  Objective:  To  find  the  largest  interval  [k,k]  such  that 
x(t)  >  0  ,  y(t)  >  0  ,  z(t)  >0  is  optimal  if  the 
machine  stock  is  between  (k,k]  . 

(2)  This  method  is  similar  to  the  dual  simplex  method.  By  maintaining 
the  dual  feasibility  and  complementary  slackness,  we  try  to  find 
the  largest  interval  such  that  the  primal  feasibility  is  also 
satisfied. 

0 

(3)  z(t)  =  k  -  r  J"  y(s)ds  t  <  0 

t 

0 

c(t)  =  l-y(t)-z(t)+q*z(t)  =  l-k  +  qk-y(t)-r(q-l)  J*  y(s)ds  . 

t 

Step  1:  Solve  equation  U'(k*q)*e^*^  with  boundary 

condition  y(0)  =  l~k  . 

Step  2:  Compute  z(t)  . 

Find  T  (T  <  0)  such  that  z(T)  =  k^  .  (T  always  exists, 
for  as  T  ^  ^  implies  the  equation  in 

Step  1  has  no  solution.) 

Step  3:  If  for  all  t  €  [T,0)  ,  0  <  y(t)  +  z(t)  <  1  ,  y(t)  >  0  ,  z(t)  >0 
Let  k  •*-  kg  ,  A  -T  .  Return. 

Otherwise,  go  to  Step  4. 

Step  4:  Find  smallest  f  (f  <  T)  ,  such  that  for  t  G  (T,0) 

0  <  y(t)  +  z(t)  <  1  ,  y(t)  >  0  ,  z(t)  >  0  . 

Let  k  *■  z(T)  ,  A  -T  .  Return. 


PROBLEM  II  (k) 


(1)  Objective:  To  find  the  largest  interval  [k,k]  such  that 
x(t)  =  0  ,  y(t)  >  0  ,  2(t)  >0  is  optimal  if  the  machine 
stock  is  between  [k,k]  . 

(2)  This  method  is  similar  to  the  simplex  method.  By  maintaining 
the  primal  feasibility  and  complementary  slackness,  we  try  to 
find  the  largest  interval  such  that  the  dual  feasibility  is 
also  satisfied. 

(3)  x(t)  -  0  ,  y(t)  =  (l-k)e"’^’'  ,  2(t)  =  1  -  (l-k)e“’'’'  , 
c(t)  =  q[l  -  (1  -  k)e  . 

s(t)  =  -U'(c(t))  +  ^rq  /  e  (c(s))ds  +  U' (q*k)J  . 

Remark;  s(t)  is  the  solution  of  the  following  integral  equation  with 

boundary  condition  s(0)  =  0  . 

00 

e"'^'^-[U'(c(t))  +  s(t)]  =  r  /  e"“'‘[(i-l)U'(c(x))  -s(x)]dx  . 

t 

Step  1:  Find  T  (T  <  0)  such  that  k^  =  1  -  (l-k)e”’^‘^  =  z(T)  . 

(T  always  exists,  for  z(T)  -*•-<»  as  T  ->■  -®  . ) 

Step  2:  If  for  all  t  €  [T,0)  ,  s(t)  >  0  . 

Let  k  k^  ,  A  •<-  -T  ,  return. 

Otherwise,  go  to  Step  3. 

Step  3;  Find  smallest  T  (T  <  T)  such  that  s(t)  >  0  for  t  G  [T,0) 

Let  k  ■*-  z(T)  ,  A  -T  .  Return. 
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ALGORITHM 

Initialize:  Find  k  decided  from  Theorem  3-2. 

Step  1:  Solve  Problem  I  (k)  . 
i  i  +  1  . 

x^(t)  x(t)  ,  y^(t)  y(t)  ,  z^(t)  z(t)  ,  k^  ic  ,  •*-  A  . 

If  k^  *  k^  ,  stop. 

Or  else,  k  k^  ,  and  go  to  Step  2. 

Step  2:  Solve  Problem  II  (k)  . 
i  i  +  1  . 

x^(t)  0  ,  y^(t)  y(t)  ,  z^(t)  z(t)  ,  k^  k  ,  A^  A  . 

If  *  ^0  ’ 

Or  else,  k  k^  ,  and  go  to  Step  1. 

To  find  the  optimal  solution,  we  have  to  reverse  the  order  of  the 
indices.  Let  N  +  1  be  the  last  index. 


old  indices -*-k. 


k. 

1 


k.  -  1 
1 


VH+hH+hH+h- 


new  indices -*-k. 


Vl-i  k^^2-i 


A, 


\n-minnn\ 


k„  N+l  1 

N 


.  2 :  i  i  N  +  i 
Vi-i .  1 ;  1 :  N  +  i 


i  „N+2-i.  .  i,  ,  N+2-i,  ,  i,  .  N+2-i,  , 

x  ■>-  X  (t)  ,  y  (t)-*-y  (t)  ,  z  (t)'^z  (t)  , 

2  <  i  <  N  +  1  . 
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Define  T. 

1 


y  A .  ,  1  <  i  <  N 
i=l  '  =  = 


Tq  =  0  .  Then 


(1)  x^(t)  ,  y^(t)  ,  z^(t)  is  optimal  for  t  £  1  $  i  ^  • 

(2)  x(t)  =  0  ,  y(t)  =  (1-^)0'*^*^  ,  2(t)  =  1  -  (1 -k)e"'^‘^  is  optimal 

for  t  >  T  . 

=  n 
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§3.  Quadratic  Utility  Function 

■) 

(M  *  X ) 

When  the  utility  function  is  quadratic,  U(x)  - - 2 - >  ^ 

we  have  a  simple  optimal  solution.  First  we  would  like  to  compute 


k  . 


=  M  -  q  +  q(l-k)e 


-rt 


s(t)  =  -h  rq  J"  e 

t 


^  r(q  -  1)  - 
a  +  r 


_  g  -  r(i  -  1)  +  r 


(M-q) 


a  +  2r 


q* (1  -  k)e 


■rt 


Case  i;  M  >  q 


If  a  -  r  (q  -  1)  r  <  0  ,  then  k  =  0  . 

If  a-r(q-l)  -f-  r  >  0  ,  let  k  satisfy  — ~  (M-q)  - 

ct  "T  r 

a - ^^*3  ~  ^  ^  q(l  -  k)  =  0  ,  then  k  =  max  {0,k}  . 

a  -I-  2r  ^ 


Case  ii;  M  =  q 

If  a-r(q-l)  -t-  r  <  0  ,  then  k  =  0  . 

If  a-r(q-l)  r  >  0  ,  then  k  =  1  .  (By  Property  2,  this  means 

we  never  reach  this  situation.) 


Remark: 

a  -  r(q-l)  -t-  r  <  0  is  not  a  sufficient  condition  for  k  =  0 
in  a  general  concave  utility  function.  If  k  =  0  ,  then  we  have 
x(t)  *  0  for  all  t  .  If  k  >  0  ,  we  would  like  to  find  the  optimal 
program  by  using  Algorithm. 


Solve  Problem  I  (k)  : 
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=  M  -  1  +  k  -  q*k  +  y(t)  +  r(q-l) 


j  y(s)ds 


t  <  0  . 


Solve  U’ =  ij/(t)  with  y(0)  =  1  -  k 


equation 


r(q  - 1) 


a  +  r 


(M-i)  -  ^  r  - 

^  a  +  2r 


q(l  -  k) 


y(t)  ,  we  get 


Using  the 
0  to  simplify 


y(t)  =^^-~.H)[a-r(q-l)  +r]^  ^  ^r(q-l)t  ^  (M-q-k)[a  -  r(H -1)1  ^fa-r(Q-l^1t 
[a  -  2r(q  -  1)  ]  (2r +a)  a-2r(q-l) 

Fact  1: 

y(t)  is  decreasing  and  y(t)  >  0  for  all  t  <  0  , 

Proof : 

a  -  r(q  -  1)  <  0  Implies  y(t)  <  0  ,  thus  y(t)  is  decreasing. 

Also  y(0)  «  0  ,  thus  y(t)  >  0  for  t  <  OB 

Fact  2: 

y(t)  +  2(t)  is  a  concave  function. 

Proof ; 

^2 

^  [y(t)  +  z(t)]  =  ^  [y(t)  +  z(t)]  <  0  (for  a-r(q-l)  <  0 
and  a-r(q-l)  +  r  >  0)  B 

Fact  3: 

y(t)  +  z(t)  <  1  for  all  t  <  0  . 


Proof ; 


y(t)  +  z(t)  is  increasing  as  t  0  .  (Otherwise,  violate 

y(t)  +  z(t)  <  1  for  y(0)  +  z(0)  =  1  .)  By  Fact  2,  y(t)  +  z(t) 
is  increasing  for  all  t  <  0  ,  thus  y(t)  +  z(t)  <  1  for  all  t  <  OB 

From  Fact  3  and  Fact  1,  Algorithm  stops. 

Summarize  the  results  as  follows: 

(1)  If  M  >  q  : 

(i)  If  k.Q  >  k  ,  then  x(t)  =  0  ,  y(t)  =  (l-kQ)e  , 

“IT  t 

z(t)  =  1  -  (l-kQ)e  is  optimal. 

(ii)  If  k^  <  k  ,  x(t)  >  0  ,  y(t)  >  0  ,  z(t)  >0  is  optimal 

if  t  <  T  ,  where  k(T)  =  ic  . 

x(t)  =  0  ,  y(t)  >  0  ,  z(t)  >0  is  optimal  if  t  >  T  . 

(2)  If  M  =  q  : 

(i)  If  a-r(q-l)  +  r  >  0  (k  =  1)  ,  then  x(t)  >  0  , 

y(t)  >  0  ,  z(t)  >  0  is  optimal  for  all  t  . 

(ii)  If  a-r(q-l)  +  r  <  0  (k  =  0)  ,  then  x(t)  =  0  , 

y(t)  =  (l-kQ)e  ,  z(t)  =  l-(l-kQ)e  is  optimal. 

Example  1: 

q=2,a*0.9,r=l,kQ=0,  and  the  utility  function  is 
(5  -  x) ^  - 

U(x)  * - 2 -  •  r(q-l)-a  =  0.1  >  0  ,  a-r(q-l)  +  r  =  0.9  >  0  . 

0.1  0.9 

X  3  -  y-g  x2x(l-k)=0  implies  k  i  0.73  .  We  can  compute 
y(t)  ,  z(t)  ,  x(t)  ,  when  k(t)  <  0.73  ,  using  the  formula  derived 
before ■ 


Let  T  be  such  that  x(T)  +  2*z(T)  =  1  ,  and  z(T)  =  k  . 


Example  2: 

Consider  a  modified  quadratic  utility  function  U(x)  =  4x  -  12 

/  c  _  V  1  ^ 

if  X  <  1  ,  U(x)  =  -  — 2  X  >  1  .  From  Example  1,  if 

k  >  k(t)  >  k  (k(t)  is  the  capital  stock  at  time  t) ,  x(t)  >  0  ; 
if  k(t)  >  k  ,  x(t)  =  0  .  If  k(t)  <  ic  ,  the  output  is  located  on 
the  linear  part  of  the  utility  function,  thus  x(t)  =  0  (by  Property  5, 
Corollary  2). 

The  following  pictures  show  the  discontinuous  behavior  of  the 
optimal  solution. 


Remark : 

We  can  remodify  Che  utility  function  to  be  strictly  concave  by 
U(x)  =  4x  -  12  -  e-(x  -  1)^  ,  e  >  0  if  x  <  1  .  As  e  >  0  ,  U"(x)  ^  0 
if  X  <  1  .  By  Corollary  1,  Property  5,  x(t)  >  0  for  all  t  ,  where 
k(t)  <  k  ,  is  not  optimal. 
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CHAPTER  IV 

GENERAL  CONCAVE  UTILITY  FUNCTION  -  MANY  MACHINES  CASE 

In  this  chapter  we  generalize  the  idea  of  Chapter  2  and  Chapter  3 
to  obtain  the  asymptotic  result.  Also,  we  give  an  example  which  shows 
the  off-beat  nature  of  this  general  problem.  Again,  we  list  the 
stationary  program  for  completeness. 

§1.  Stationary  Program 
Theorem  4-1: 

r.q. 
i  '  i 

1  >  max  — ; -  is  the  sufficient  and  necessary  condition 

=  l<i<n  “  '^i 

for  the  existence  of  an  optimal  stationary  program® 

r.q. 

1  ^  i 

From  now  on,  assume  1  <  max  -  . 

l<i<n  1 

§2.  Nonstationary  Program 

To  begin  with,  we  prove  some  simple  properties. 

Property  1: 

w(t)  is  continuous. 

Proof : 

We  only  need  to  consider  those  T  >  0  such  that  there  exist 
i  ,  j  (i  j)  and  y^^Ct)  >  0  for  t  €  (T-e,T)  and  y^  (t)  >  0 
for  t  G  (T,T  +  e)  . 
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w(T) 

w(t) 


00 


If 

t  E 

(T  -  e,T) 

w(t)  =  J 

w^ (s)ds 

00 

A 

If 

t  e 

(T,T  +  e) 

w(t)  =  r.  J 

w.  (s)ds 
3 

J  t 

00 

A 

If 

w(T) 

>  r .  1  w. 
^T  ^ 

(s)ds  ,  then 

as  t  -► 

=  r.  I  w.(s)ds  ,  contradiction.  Thus 
3  3 

is  continuous  ■ 


oo 

,  w(t)  >  r.  f  w. (s)ds 

*  1  •r  i 

t 

00 

T  ,  w(t)  r .  J*  w .  (s)d 
1^1 

OO 

w(T)  ~  f  w^(s)ds  and 


s  < 


Remark; 

This  property  gives  us  the  boundary  condition. 

Property  2: 

1^1  (t)  is  continuous. 

Proof ; 

If  we  can  prove  that  for  any  sequence  t^  T  ,  such  that  ^ 

or  there  exists  i  such  that  z^(t^)  ->■  6  and  ^  s  <  i-1  , 

we  have  >i'(T)  ,  then  i|)(t)  is  continuous. 

The  proof  is  similar  to  that  of  Property  4,  Chapter  3  I 

Corollary : 

If  u(x)  is  strictly  concave,  then  (1)  c(t)  is  continuous; 

(2)  x(t)  and  z^(t)  are  ail  continuous  1  <  i  <  nB 

Below  we  study  the  asymptotic  behavior.  We  follow  the  notation 

i  N 

of  Chapter  2,  and  let  L  =  ^  ^  ^  ® 
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For  any  best  machine  j  ,  k  =  (k^^,  ....  k^)  g  L  ,  define  the  following: 


:(t,k)  =  k..q.  +  q.^1  -  k.j.(l  -  e  ^  ) 


J-(a+r.)s  (a+r.)t 

e  ^  U' (c(s,k) )ds‘e  ^ 


G  =<K 


K  G  L  ,  and  r 


/-(a+r.)s 

e  ^  -U'  (c (s ,k) )ds' 


(a+r^ ) t 


qf^^l-U'(c(t,k))  >  0 


oo 

J 


e  *0'  (c  (s  ,k)  )dsT 


r .  t 


/-(a+r.)s 

e  ^  *11' (c(s,k))ds' 


e  ^  •q,(r.-r.)  >  0  for  i  £  I  -  {j} 
J  J  1  = 


for  all  t  >  0  . 


Given  any  initial  stocks  k^  of  machine  i,  1  <  i  <  n  ,  let 
k  =  (k^,  kjj)  ,  and  we  have  the  following  theorem. 


Theorem  4-2: 


If  KEG,  then  the  best  program  with  respect  to  j  is  optimal, 


If  N  =  n  ,  then  replace  by  1  ;  *  implies  s(t,k)  >  0  for  all  t 


The  output  rate  function  is  c(t,k)  .  Define: 


w(t)  =  e“°‘'^-[U*(c(t.k))  +  s(t,k)] 

„^(t)  =  e~“^-[(q.  -  l)U'(c(t,k))  -  s(t,k)]_j^  1  ^  i  • 

If  i  e  I  , 

(q^  -  1)U’ (c(t,k))  -  s(t,k) 

J-(a+r.)s  _ 

e  ^  -U' (c(s,k))ds-e 

t 

>  q  .u’(c(t,k))  -  u’(c(t,k))  (U’(c(t,k))  is  decreasing) 

>  0  (definition  of  I) 

implies  w^(t)  =  e  • [ (q^  -  1)U' (c(t,k) )  -  s(t,k)]  . 

If  i  €  1  , 

(q.  -  1)U' (c(t,k))  -  s(t,k) 

1 

/-(a+r.)s  _  (a+r  )t 

e  ^  -U’ (c(s,k))ds-e 

t 

<  0  (definition  of  G) 

implies  w^(t)  *  0  • 

00 

w(t)  “  r  f  w  (s)ds  (by  the  result  in  the  one  machine  case). 


If  i  ^  I  , 


w(t) 


w^(s)ds  =  w(t) 


>  0  , 


If  i  €  I  , 


w(t) 


00 

-  r.  I  w. 


(s)ds 


00 

/ 


e  (c  (s,k)  )dsT 


oo 


- (a+r . )s 


•U ' (c (s ,k) )ds’q . (r  ,  - 


:  J 


>  0  (definition  of  G)  . 


Complementary  slackness  is  also  satisfied.  Thus  the  best  program 
with  respect  to  j  is  optimal® 


Let  a 


N  /  ^  \ 

■j,  j,  hj- 


then  c(t,k)  = 


-r .  t 
3 


a  -  b*e  ,  U'(a)  <  U'(c(t,k))  <  U'(a-b)  . 


Corollary  1: 


If 


/ 


U' (a  -  b) 


q.(r.  -  r.)‘a 
J _ L--  ^ 


r  .q. 

U'(a)  *7 — r^— 4 -  ,  U'(a)  ,  . 


if  i  £  I  -  tj)  and  r^  >  r^  , 


<min/  ^ 

I  (a  +  r  )r  -(q  -  q  ) 

/u'(a) - H — - — r* - ^  if  i^I 

I  ''  n  (r  —  ri^r* 


{j}  and  r 


then  the  best  program  with  respect  to  j  is  optimal. 


Proof ; 


r.q.  I  e 

rj  J 

t 


- (a+r . ) s 


J  .TT' 


(ot+r . )  t 


U' (c(s,k))ds*e 


-  'Jn+i'U'  (c(t,k)) 


^  ^  U'(a)  -  (a  “  b)  >  0  (by  assumption) 


If  r.  >  r.  and  i  6  I  -  {j} 
J  1 


00 

/ 


e  (c(s,k)  )ds*r 


/-(a+r.)s  r.t 

e  ^  ‘U' (c(s,k) )ds*e  ^  ‘q^'Cr 


ru'(a-b)*r.»(q.  -  q.)  U'(a)*q.*(r.  -  r.)"] 

>  - i - ] - h-  +  - — J - 3 - 

=  L  “  «  +  r .  J 


e  (q^  >  q^) 


>  0  (by  assumption) 


If  r 


^  >  r J  and  i  G  I  -  { j } 


00 

/ 


e  (c(s,k))ds 


/- (a+r . ) s  _  r.t 

e  ^  ’U' (c (s ,k) )ds* e  ^  •q^-(r 


j"!!'  (a)  -r^*  (q^  ”  "^i^  U’  (a  -  b)*q^*  (r^  ~  ’^i^  1 
=  L  ®  ^  a  +  r.  J 


-at  /  , 

e  (q^  >  q^) 


>  0  . 


Thus  K  G  G  and  the  best  program  is  optimal  I 


Corollary  2: 

If  U'(a)  >  0  and  — ^  for  i  j  ,  then  when 

M  ^  ^ 

I  k,  ^  1  ,  ic  G  G  . 
i-1  ^ 
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Proof : 


r . q .  r .q  . 

ct  +  ^N+1  (a  + 

:  i  G  I  -  {j}  and  ’^j  ^  ’ 

q.(r.  -  r.)a  -  (a  +  r.)r.(q.  -  q.)  =  r.q.(ot  +  r.)  -  r.q.(a  +  r.) 
J  J  1  j  1  1  i^i  j 

(r.q.  r.q.  \ 

a  r.  a  +  r J  ’ 


>  0 


q .  (r  .  -  r .  )a 

implies  7 — T*^ — ^ - y-^ - r  >  1  . 

(a  +  r  Jr^(q^  -  q^) 


If  i  €  I  -  {j}  and  » 


(a  +  rj)r^(qj  -  q.)  -  q^(r^  -  r^)a  =  ^  °  ’ 

(a  +  r,)r  • (q.  -  q.) 


implies 


q.(r.  -  r.)a 
J  1  1 


N 


>  1  . 


N 


As  \  k  -*•  1  ,  U'(a)  -*■  U'(a-b)  ;  thus  as  \  k.  -*■  1  ,  Corollary  1 
i=l  ^  i=l  ^ 

is  true,  which  implies  K  G  G  ■ 


Remark : 

If  the  best  machine  is  not  unique,  we  can  use  Taylor  Expansion 
to  get  the  similar  result. 


Lemma  1 : 

As  t  ,  c(t)  -*■  constant,  2^(t)  -*•  constant,  yj(t)  -►  0  , 


x(t)  -♦  constant. 
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Proof : 

Since  k^(t)  is  increasing  and  bounded,  k^(t)  "*■  as  t  , 

n 

which  implies  y.(t)  -*■  0  ,  ^  z.(t)  +  x(t)  ->■1  as  t  . 

^  i=l  ^ 

Since  q,  >  q„  >  . . .  >  q  >  1  ,  so  x(t)  and  z.(t)  ->■  constant 
1  z  n  1 

n 

as  t  .  c(t)  =  x(t)  +  q.*z.(t)  -*•  constant  as  t  ->■ 

i=l  ^  ^ 


Theorem  4-3: 

N 

As  t^“>,  5!  k.(c)-*-l. 

i=L  ^ 


Proof : 


For  any  e  >  0  ,  there  exist  T,P, ,P.,  P  ,  P  >  0  ,  such  that 

1  /.  n 

n 

x(t)  >  P  ,  z.(t)  >  P.  1  <  i  <  n  and  J  P.  +P>l-e  for  all 
^  ^  i=l  ^ 

t  >  T  .  After  time  T  we  can  solve  the  original  problem  in  the 
following  two  steps: 


Step  1: 
Step  2: 


Put  P^(P)  units  of  labor  to  work  with  machine  i 


Consider  a  new  problem  with  1 


-(I 


units 


(alone) . 
of 


labor:  let  c(t)  be  the  output  rate  function,  U(x) 
be  the  utility  function,  where  U(c(t))  = 

ujico  .  p,,,  *  p)  . 


As  e  -►  0  ,  u(c(t))  -►  constant  for  all  t  ,  thus  the  optimal 
program  for  Step  2  is  the  best  program  with  respect  to  j  ,  which 
implies  that  we  only  build  machine  j  as  t  .  Combining  Step  1 
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and  Step  2,  we  find  that  as  t  ->•  “  we  only  build  machine  j  .  From 

the  proof  of  Theorem  4-2,  we  must  have  z^(t)  =  0  i^I  as  t-*-®, 

N 

which  implies  J]  k  (t)  -►1  as  t  ■ 
i=l  ^ 

Corollary: 

If  U'(a, )  >0  ,  ^  ^  ^  ^  for  i  /  j  ,  then  for  t 

sufficiently  large,  the  best  program  with  respect  to  j  is  optimal ■ 

Property  3: 

i(/(t)  is  decreasing. 

Proof : 

By  contradiction.  If  not,  let  (T-e,T)  be  the  last  time  interval 
such  that  iKt)  is  increasing.  (T  exists,  for  4'(t)  is  continuous 
and  as  t  -►  “>  the  best  program  is  optimal.) 

There  are  two  possibilities:  x(t)  >0  or  0  <  z^(t)  <  k^(t) 
for  some  i  .  If  x(t)  >  0  and  y.(t)  >  0  ,  t  €  (T-e,T)  , 


-at 


w(t)  =  e  -iJiCt)  =  r^^  I  Wjj^(s)ds 


00 

/ 


] 


J 


<  I  Wjj^(s)ds  -(•  -  1) ‘iijCT)  |  e  “^ds 


(iKt)  is  decreasing  for  t  >  T  and  < 


-at 


'i'(t)  •  (q  -  1))  . 


implies  a  <  r^(q^-l)  . 


As  t  ->■  T_  ,  -  i*-(T)-e"“‘’^ 

But  4)(t)  =  c*e  (by  Property  5,  Chapter  3)  implies 

is  decreasing,  contradiction.  In  the  case  of  0  <  2^(t)  <  k,^(t)  , 
the  proof  is  similar  ■ 

Corollary  1: 

If  y^^Ct)  >  0  then  a  <  * 

Proof : 


e  <  w(t)  ■ 


Corollary  2: 

If  u(x)  is  strictly  concave,  then  c(t)  is  strictly  increasing® 


§3.  Example 
Example  1: 

The  following  two-machines  problem  shows  that  we  build  different 
machines.  (Time  [0,T)  we  build  machine  2,  time  [T,“>)  we  build 
machine  1.) 

2 

U(x)  =  -  ,  a  *  0.9  ,  *  1  ,  “  50  ,  r^  “  4  ,  ,  q  * 


"l^l 
a  +  r. 


.  _!2!2_ 
'*2  a  +  r. 


implies  machine  1  is  best,  machine  2  is  good. 


1  -  -  .  -2 

And  we  can  find  the  exact  dual  price,  using  the  technique  in  Chapter  3, 


such  that  there  exist  T 


,  k  >  0  (v 


1  -  e 


S  0.09321  and 


-r  (t-T)  -r  (t-T) 

x(t)  =  0  ,  yj^(t)  =  (l-k)e  ,  Zj^(t)  =  1  -  k  -  (l-k)e  , 

Z2(t)  =  k  ,  y2(t)  =0  is  optimal  for  t  >  T  ;  x(t)  «  0  ,  y^^Ct)  =  0  , 

-r2t  -r2t 

Zj^(t)  =  0  ,  y2(t)  “  e  ,  Z2(t)  =  1  -  e  is  optimal  for  t  <  T  . 


Example  2: 

This  two  machines  example  shows  that  we  build  machine  2  in  some 


time  interval  (Tj^,T2)  ,  but  as  t  -*•  “>  we  never  use  machine  2. 
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M  =  51  ,  r^  =  1  ,  =  50  ,  k^  =  k,  =  0  ,  a  =  0.9  ,  q,  =  -j— g  , 


r2  =  some  large  number  (decided  later) ,  u(x)  =  - 


(51  -x)- 


1  , 


’^l‘’l  50  40 


^2^2 


■  t:?  "  TT?  ■  ’2  2  '2<'>  -  0  ' 


(machine  2  is  not  good),  a  -  rj^(q^  -  1)  +  r^  =  0. 9  -  49  +  1  <  0  , 


-rit  -r^t 

so  if  we  can  prove  that  x(t)  *  0  ,  y^^Ct)  =  e  ,  Zj^(t)  »  1  -  e 


y2(t)  =  0  ,  22^*-^  *  ®  optimal,  then  there  must  exist  some  time 


-r,t 


interval  in  which  we  build  machine  2.  If  x(t)  =  0  ,  yj^(t)  =  e  ^  , 


"^1^  -at 

2^(t)  =  1  -  e  is  optimal,  then  W2(t)  =  e 


/  "l"!  \ 

-I- ^->2  -  ■ 


('’2  - 

»  '  /I 


If 


'’2 


a  +  r,  ’  ^*2  a  +  2r, 


>  0 


and 


then  there  exists 


(“12  -  ''l(‘'2  -  “  • 

T>0  such  that  ^92  “  ®  *'’l’ ^‘^2  “  °  > 

which  implies  W2(t)  *0  t  >  T  ,  W2(t)  >0  t  <  T  . 


r2^  W2(s)ds  ^  ^2  J  *  r2*k  where  k  »J  W2(s)ds  >  0 


T 

2 


1 

2 


T 

2 


r„ 


as 


-►  00 


Contradiction  to 


q.  -  -i  —  =  -5  <  0 

^2  a  + 


c 

^2  2r  +  a  ~ 


Check: 


20  -  16.3  a:  3.7  >0  . 
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CHAPTER  V 

VARIANT  OF  THE  MODEL 


In  this  chapter  we  introduce  the  depreciation  factor  6  into  the 
model.  Difficulty  arises  here  because  the  capital  stocks  are  no  longer 
monotonic  over  time  in  the  many  machines  case.  For  simplicity,  we 
concentrate  on  the  one  machine  case. 

The  model  with  the  depreciation  factor  5  is: 


00 

/ 


-at 


<P)  max  I  e  •U(x(t))dt  s.t.  A*x(t) 

*0 

x(t)  >  0  . 


(s)ds  +  b(t)  , 


00 

I 


CD 

/-6(s-t) 
e  ‘w 


<D)  min  I  w(t)*b(t)dt  s.t.  w(t)A  > 

*D  t 

e  °‘^*VU(x(t))dt  ,  w(t)  >  0  . 


(s) ‘Eds  + 


In  the  one  machine  case,  if  we  compare  the  depreciation  and  non- 
depreciation  cases,  the  differences  are  the  following  two  equations: 


I 


z(t)  5  r  I  e  ®^*y(s)ds  +  k^^e 


w(t)  >  r 


00 

f  -6(s-t) 

J  *''l 


(s)ds 


The  following  program  is  called  contraction:  x(t) 
y(t)  »  0  ,  z(t)  *  kQ-e"*^^  . 


l-kQe 


-6t 
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1 


Theorem  5-1 : 


1  >  — ; — ^ — ; — 7  is  a  necessary  condition  for  the  existence  of  an 
=  a  +  r  +  6 

optimal  contraction  program. 


Proof ; 


*  5 1  “  5 1 

c(t)  =  q*kQ*e  +  1  -  kQ*e  ,  i|)(t)  is  increasing. 
Define:  w(t)  =  e  “*'-4)(t)  ,  =  e  (q  -  l)i|; (t)  . 


00 

<  w(t)  -  r^  e 


<  e  '^*'*ij/(t)  -  r  (q  -  l)ijj(t) I  e 


QO 


(4)(t)  is  increasing) 


01  T  0 


implies  a  +  6  -  r(q  -  1)  >  0  ,  1  >  ^  ^ 


Remark : 


(1)  (*)  is  the  sufficient  and  necessary  condition  for  the 
existence  of  an  optimal  contraction  program. 

(2)  In  the  linear  utility  function  case,  the  above  condition 
is  also  sufficient. 


From  now  on,  assume  1  <  — — ^ - — r  and  U'(q)  >  0  . 

a  +  r  +  6  ^ 

In  the  following,  we  prove  an  asymptotic  result  similar  to  the 


nondepreciation  case. 


Lemma  1 : 


If  x(t)  =  0  ,  z(t)  =  r  /  e  ®^*y(s)ds  +  k*e  ^  ,  x(t)  +  y(t)  + 

0 

6  /  r  -(r+6)t  r  _ 

z(t)  =  1  ,  then  y(t)  =  j.  ^  ^  "*■  (r  +  6  “  r  +  6 


-(r+6)t 


Proof : 


r  f  e  ®^y(s)ds  +  y(t)  = 


,  ,  “5t  .  ,  . 

1  -  k-e  ,  implies 


Define: 


k)  =  -U’(c(t,k))  +  rq  J  (c(x,k))dx*e^“'''’^'^'^^’^ 

t 

k)  =  U’(c(t,k))  -  vj  (c(x,k))dx*e''“'''’^''’'^^’^ 

t 

(k  I  0  <  k  <  1  .  s(t,k)  >  0  ,  9(t,k)  >  0  for  all  t  >  0} 


a  +  r  +  6 


a  +  r  +  5 
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Lemma  2 : 

s(t,k)  is  the  solution  of  the  following  equations: 
w(t)  =  e  “*'*[U'(c(t,k))+s(t,k)]  ,  w^(t)=e  “*'*[(q-l)U'(c(t,k))-s(t,k)], 

oo 

w(t)  =  r  J*  e  ^  ^^•w^(s)ds  . 
t 

Proof : 

It  is  a  straightforward  calculation  ■ 

Theorem  5-2: 

If  k  €  G  ,  then  x(t)  =  0  ,  y(t)  =  |  , 

=  7^  "  (tTT  ■  is  optimal. 

Proof : 

Define : 

{w(t)  =  e  [U' (c(t,k) )  +  s(t,k)] 
w^(t)  =  e  “^'[(q -l)U'(c(t,k))  -  s(t,k)]_|_  , 

00  - 

iTi  I  -(a+r+6)s  ,  ,  -  .  ,  (r+a+6)t 

U  (c(t,k))  -  r  I  e  *11  (c(s,k))ds*e 

J  + 

=  e  ““^^'q-eCt.k)  =  e  (  (q  -  1)U  '  (c  (t ,  k) )  -  s(t,k)]  . 

By  Lemma  2  and  the  definition  of  G  ,  the  rest  of  the  proof  is 
straightforward  I 

Corollary : 

If  k-,  -►  — 7— r  ,  then  k.  G  G  . 
u  r  +  0  U 


W^(t) 


-at  - 

e  •  q* 
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Proof ; 


As  k 


0  r  + 
implies  €  G  ■ 


,  s(t,kQ) 


>  0  ,  eCc.kp) 


Lemma  3 ; 


1  I  “<5t  ,  r  -<5(t-s)  /  V  .  .  ^  . 

k(t)  =  k-*e  +  r  J  e  •y(s)ds  is  monotonic. 


Proof : 

If  k(t)  is  not  monotonic,  say  increasing  in  (T-e,T)  ,  decreasing 
in  (T,T  +  e)  .  Since  k(t)  is  continuous,  some  moment  right  before  T 
should  have  the  same  capital  stock  as  some  moment  right  after  T  ,  thus 
the  same  behavior.  (See  the  illustration  below.) 


Contradict  to  the  fact  that  k(t)  is  increasing  before  T  and  k(t) 
is  decreasing  after  T  ■ 

Lemma  4 : 

As  t  ->■  ®  ,  c(t)  -♦  constant,  ■*  constant. 


Proof : 

As  t  -*■  ®  z(t)  =  k(t)  ■*  constant  (k(t)  is  monotonic  and  bounded); 
y(t)  -*■  constant  (k(t)  -*•  constant); 


x(t)  -*■  constant  (x(t)  +  y(t)  +  z(t)  =  1)  ; 

c(t)  =  x(t)  +  q*2(t)  -*•  constant,  liiCt)  -*■  constant  ■ 

Lemma  5 : 

If  x(t)>0,y(t)>0,z(t)>0,  then  ip(t)  =  c*e^'^"'’'^ 

Proof : 

Straightforward  ■ 

Lemma  6 : 

As  t-*-”,y(t)-*-6>0. 

Proof : 

By  contradiction.  If  y(t)  -*•0  as  t  ,  then  z(t)  0  , 
x(t)  -*•  1  as  t  .  If  y(t)  #0  as  t  ,  contradiction  to 
Lemma  4  and  Lemma  3.  If  y(t)  =  0  as  t  ,  contradiction  Co 
Theorem  5-1 B 

From  Lemmas  4,  5  and  6,  we  know  for  t  sufficiently  large, 

x(c)  =  0  .  Thus  k(t)  -*■  — 7—7  ,  k(t)  EG  as  t  .  We  state  the 

r  +  0 

above  results  in  the  following  theorem. 

Theorem  5-3: 

(1)  If  k„  <  — 7— r  Chen  k(t)  is  increasing,  and  for  t 

U  r  +  0 

sufficiently  large  k(t)  ->■  ,  x(t)  =  0  . 

(2)  If  k_  >  — 7—7  then  k(t)  is  decreasing,  and  for  t 

sufficiently  large  k(t)  -*■  ,  x(t)  *  0  ■ 
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